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ABSTRACT 

We present timc-resolved spectroscopic and polarimetric observations of the AM Her system EU 
Cnc. EU Cnc is located near the core of the old open cluster Messier 67; new proper motion mea- 
surements indicate that EU Cnc is indeed a member of the star cluster, this system therefore is useful 
to constrain the formation and evolution of magnetic cataclysmic variables. The spectra exhibit 
two-component emission features with independent radial velocity variations as well as time-variable 
cyclotron emission indicating a magnetic field strength of 41 MG. The period of the radial velocity 
and cyclotron hump variations are consistent with the previously-known photometric period, and the 
spectroscopic flux variations are consistent in amplitude with previous photometric amplitude mea- 
surements. The secondary star is also detected in the spectrum. We also present polarimetric imaging 
measurements of EU Cnc that show a clear detection of polarization, and the degree of polarization 
drops below our detection threshold at phases when the cyclotron emission features are fading or not 
evident. The combined data are all consistent with the interpretation that EU Cnc is a low-state 
polar in the cluster Messier 67. The mass function of the system gives an estimate of the accretor 
mass of AfwD > 0.68Mq with Mwd ~ O.SSMq for an average inclination. We are thus able to place 
a lower limit on the progenitor mass of the accreting WD of > 1.43Mq. 

Subject headings: white dwarfs - novae, cataclysmic variables - Stars: individual; EU Cnc - open 
clusters and associations: individual: Messier 67 - Accretion, accretion disks 



1. INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary sys- 
tems in which a white dwarf (WD) is accreting material 
from a low-mass companion star. If the WD has a suffi- 
ciently strong magnetic field, the formation of an accre- 
tion disk is inhibited, and material accretes directly onto 
one or more magnetic poles of the white dwarf. These bi- 
naries are known as AM Her systems or polars, after the 
high fraction of polarized light detected in the systems. 

Ikurtis.williams@tamuc.edul 



The identification of CVs in star clusters provides in- 
teresting constraints on the formation and evolution of 
the interacting system. As all members of an cluster are 
coeval, the total age of the system is known, as is the 
system's distance and metallicity. Further, if the mass 
and effective temperature of the white dwarf can be de- 
termined, then a limit on the white dwarf's progenitor 
mass can be derived via the same metho ds used to con- 
struc t the initial-final mass relation (e.g.. lWilliams et al.l 
l2009f ). Because the white dwarf is likely re- heated to 
some extent by the ongoing accretion, the constraint on 
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the progenitor niass would be strictly a lower limit. 

CVs are quite common in globular clusters (e.g. 
IMareon et all [ToSlt IGrindlav et all ITOOl . The glob- 
ular CV p opulation tends to b e cen trally concen- 
trated (e.g. IGrindlav et all [TOOl [20Q1 . This resuh 
is likely explained by the formation of tight binaries 
by stellar encounters in globular cluster cores (e.g. 
iPodsiadlowski et al.ll200l iPoolev et al.ll2003D . As such, 
globular cluster CVs are excellent tracers of a globular's 
dynamic history, but these CVs may not shed much light 
on the formation and evolution of CVs in the galactic 
field where stellar encounters are rare. 

Open star clusters may therefore be a more useful lab- 
oratory for studying CV evolution. The stellar densities 
are far lower than in globular clusters, even in the cluster 
core, and dynamical simulations suggest that CV forma- 
tion is not enhanced by the stellar encounters that do 
occur dSh ara fc Hurlevlf200l . though a small fraction of 
open cluster CVs may s till be formed by stellar exchanges 
(jShara fc HurlevI [20061 ) . Open clusters also span a wide 
range of ages, metallicities, and stellar masses, raising the 
potential to study how these parameters impact CV for- 
mation and evolution more precisely than possible from 
studies of field CVs. 

Unfortunately, the number of CVs in open clusters is 
small, and none are well-studied. In the ancient, dense, 
metal-rich open cluster NGC 6791, two spectroscopically- 
confirmed CVs are known ( K aluz nv et al. 199^; 
iMochejska et all 120031 ): Ide Marchi et all (2007) identify 
a suspected third cluster CV based on photometric 
properties and conclude that all three CVs are likely 
cluster members. One CV, EU Cnc, is known in the 
open cluster M67 and described in detail below. These 
four objects are the only confirmed cluster member 
CVs, and due to the large distances of the clusters 
[{m - M)v = 13.4 for NGC 6791 and (m - M)v = 9.97 
for M67], spectroscopic studies of these CVs with the 
same precision and the same techniques (such as tomog- 
raphy) as current field CV studies require significant 
time on 8 m-class telescopes and larger. 

Other candidate open cluster CVs have been sug- 
gested: .Mocheiska et al.l (j2004 12006) identify a CV in 
the field of the ~ 2 - 3 Gyr-old open cluster NGC 2158, 
though it may lie foreground to the cluster. One CV 
is identified photometrically in the field of the 3.5 Gyr- 
old cluster NGC 6253, b ut no membership information is 
available (jde Marchi et a l. 2010). The rich open cluster 
M37 (age ~ 550 M yr) has two CV c andid ates identified 
photometrically by iHartman et al.l (|2008f) . Finally, an 
X-ray source in the field of the cluster NGC 6819 (age 
~ 2 — 2.4 Gyr) has properties consistent with CVs, but 
its true nature and cluster membership are unconfirmed 
(jGosnell et al.l [201^ . T he Hyades con tains at least one 
pre-CV, V471 Tau (e.g lVauclai3ll972| ). but no actively 
accreting systems. 

1.1. EU Cancri 

EU Cnc was detected as variable object in the field of 
the old open star cluster Messier 67 by Gillilan d et al.l 
|K)91), who identified it as a likely AM Her system 
based on the similarities of its light curve to that of 
VV Pup. They determined EU Cnc has a photomet- 
ric period of 2.091 ± 0.002 hr with variations through a 
CuS04 [U band) fiher of 0.6 mag. .BcUoni et al.. (,1993.1 



detected an X-ray source coincident with the optically- 
variable source; its very soft X-ray hardness ratio is typ- 
ical for AM Her systems in the ROSAT band s. The AM 
Her n ature of EU Cnc was confirmed by Pa squini et al.l 
(|1994i). who obtained three 75-minutc optical spectra of 
the source. These spectra exhibit cyclotron humps and 
radial- velocity variable emission lines of H, He I, He II, 
and Fe II. Under the assurnption that EU Cnc is in the 
star cluster, iPasauini et al.l |l994[) conclude that the ab- 
solute optical magnitude and X-ray luminosity of EU Cnc 
are typical for low-state AM Her systems. 

Subsequent X-ray studies bv lBelloni et al.l ([1998') using 
ROSAT detected 100% modulation of the soft X-ray flux 
with a period equal to the optical period, indicative of 
accretio n onto a single magnetic po le. Chandra observa- 
tions bv Ivan den Berg et al.l (|2004D detected hard X-ray 
emission from EU Cnc, again typical for AM Her systems 
and likely due to shocks in the accretion fiow.^ 

Mor e recent time-series photometry from iNair et all 
(l2005l) again detected high-amplitude optical modula- 
tion, with magnitudes var ying from 21.6 to 20. 3 mag 
at the same 2.09 hr period oflCillil and et al.l (|1991D . The 
variation amplitude was about 30% larger than in 1991, 
with that difference likely due to different filters used in 
the two studies, V containing a large cyclotron modula- 
tion. 

Based on the body of work on EU Cnc, INair et all 
([2005) point out an interesting conundrum. The high- 
amplitude optical variability is typical for AM Her sys- 
tems in a high accretion state, while the X-ray luminosity 
of EU Cnc, under the assumption it is a member of M67, 
is typical of magnetic CVs in a low accretion state. 

In this paper, we present time-resolved spectroscopy 
of EU Cnc obtained serendipitously with the 10-m Keck 
telescope, as well as the first polarimetric measurements 
of this system. After discussing the observed phe- 
nomenology, we revisit the issue of EU Cnc's cluster 
membership. We will show that the optical variability 
is almost entirely due to cyclotron emission changing 
throughout the orbit and that the optical spectrum is 
typical of a low accretion state polar. We detect the sec- 
ondary star in the spectrum and, along with the other 
evidence, we show that EU Cnc is a member of the M67 
open cluster. 

2. TIME-RESOLVED SPECTROSCOPY 
2.1. Observations and Data Reduction 

We targeted EU Cnc serendipitously as part of a pro- 
gram to obtain high signal-to-noise spectroscopy of WDs 
in Messier 67. We obtained observations on UT 2007 Jan 
19 with the low-resolution imaging spectrometer (LRIS) 
on th e Keck I telescope (|Oke et al.l 119951 : iSteidel et al.l 
I2004D . We obtained simultaneous spectra with both blue 
and red sides of the spectrograph through a multiobject 
slit mask with slitlet widths of I'.'O. On the blue side, 
we used the 400 lines mm~^ grism blazed at 3400A for 
a resulting spectral resolution of ~ 7A FWHM. On the 
red side of the spectrograph, we used the D560 dichroic 
and the 600 grooves mm"^ grating blazed at 7500 A, for 
a resulting spectral resolution of 4.8A FWHM. The to- 
tal spectral coverage ranged from 7400A blueward to the 
UV atmospheric cutoff. 

The weather was nearly photometric through the en- 
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Figure 1. Time-ordered spectroscopy of EU Cnc. The left panels contain the entire spectrum; the gap at A 5650A is the gap between 
the blue and red arms of the spectrograph. The right panels contain close-ups of the A4686 He II line (left) and H/3 line (right); the radial 
velocity variations and variable line asymmetries are clearly visible. Phase </< = indicates the negative zero crossing of the narrow-line 
component of the emission lines. 



tire observation, and seeing was moderate at 0'.'9 FWHM. 
Seven exposures, each of twenty- minute integration, were 
taken over a ~ 2.5 hr period, with two short breaks 
for mask re-alignment. During the final exposure, the 
flux dropped dramatically for most stars, likely indicat- 
ing that the mask was slightly misaligned. 

These observations were taken prior to the installation 
of the LRIS atmospheric dispersion corrector, and due 
to constraints in slit mask design, the slitlets were ori- 
ented nearly perpendicular to the parallactic angle. We 



attempted to minimize the effects of atmospheric dis- 
persion by using a blue filter in the guider camera, but 
diminution of the UV light is severe at higher airmasses. 

We reduced the data using the onedspec and twodspec 
packages in IRAF. Overscan regions were used to de- 
termine and remove the amplifier bias. Flat-fielding on 
the blue side data was accomplished usin g a piecewise- 
smooth response function as described in Will iams et al.l 
(|2009l ) to eliminate ringing due to a sharp inflection point 
in the flat field at ~ 4200A. Cosmic rays were removed 
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from each two-dimensional spectrum using the L.A. Cos- 
mic L aplacian cosmic ray rejection routine (|van DokkumI 
120011) . Wavelength solutions on the blue side were de- 
rived from Hg, Cd, and Zn arclamp spectra; on the red 
side, Ne and Ar arclamp spectra were used. These cal- 
ibrations were obtained prior to the final mask realign- 
ment. 

Relative fiux calibration was obtained using l"-wide 
long-slit spectroscopy of the spectrophotometric stan- 
dards G191-B2B and G138-31 taken at parallactic angle 
near the start of the night. No attempt was made to 
obtain absolute spectrophotometric calibrations. 

2.2. Spectral Phenomenology 

The time-series spectra are shown in Figure [1] Qual- 
itatively, the spectra appear fairly typical for AM Her- 
type systems in a lo w accretion state, similar to those ob- 
served for HU Aqr (IGlenn et al.lll994 ) . Cyclotron humps 
are visible and variable in strength. Emission lines of H 
and He are observed and have at least two components, a 
narrow (unresolved) component and a broad component. 
The emission lines have variable radial velocities, and the 
radial velocities of the two line components are not in 
phase. This type of emission line compone nt behavior is 
often seen in polars, for example VV Pup (jMason et al.l 
I2008D whereby the narrow line compo nent phases with 
the motion of the secondary star (iMason et al.l 120081 : 
IHowell et aLll2008D . The time scales for both the cy- 
clotron hump variations and emission line variations ap- 
pear to be consistent with the known photometric period, 
though additional time-series spectra would be required 
to prove these are identical. We now quantify these phe- 
nomena. 

We estimate the strength of the magnetic field us- 
ing the spacing of the cyclotron humps in the optical 
spectra. The locations of the peaks are estimated by 
fitting a 12th-order polynomial to the continuum, ex- 
cludi ng the obvious emission lines . We use equation 54 
from iWickramasinghe fc Ferrariol (|2000( ) to get a mag- 
netic field strength of B « 41 MG. We do not see Zee- 
man split hydrogen absorption lines in the spectrum cor- 
responding to this magnetic field strength, but later we 
will see that they are likely filled in by the red continuum 
contribution of the secondary star. 

2.2.1. Radtal Velocity 

Because these data were not obtained with the goal 
of obtaining precision radial velocities, and because the 
flexure in the spectrograph was significant over the course 
of the observations, it was not possible to obtain abso- 
lute radial velocities. Attempts to use night sky emission 
lines as velocity zero points failed due to the lack of night 
sky lines in data from the blue arm of the spectrograph, 
where the majority of the CV's emission lines were lo- 
cated. 

Instead, relative radial velocities were obtained as fol- 
lows. First, the continuum was fit with a high-order 
polynomial and removed from each spectrum. The ra- 
dial velocity of the narrow component of the emission 
lines was determined by autocorrelation of each spectrum 
with that of an observation with an exposure midpoint 
close to zero phase. This determination involved some 
iterative bootstrapping with the radial velocity fitting 




Phase 

Figure 2. Radial velocity curve for EU Cnc (top panel) and a 
G-type star in a neighboring slitlet (bottom panel). Filled circles 
(blue in the online color version)indicate relative radial velocities 
from emission lines on the blue arm of the spectra; stars (red in the 
online version) indicate data from the red arm of the spectrograph. 
The solid curve (magenta in the online version) is the best-fitting 
sine wave. 



described further below. The resulting radial velocities 
(relative to the observations closest to zero phase) are 
given in Table [1] and shown in Figure [2j The radial ve- 
locities at « 0.59 are likely measuring different motion 
than the other points, as the narrow component of the 
emission lines disappears during these observations. 

As a cross-check, we used the same method to deter- 
mine the radial velocity of a G-type star targeted in a 
neighboring slitlet. The radial velocity of this star should 
be constant over the set of observations, and although a 
slight positive velocity offset may exist, it is of low ampli- 
tude compared with the variations observed in the CV. 

The phase and amplitude of the radial velocity varia- 
tions were determined by fitting a function of the form 

Vrei — if sin[27r((/) — 00 + 0-5)] + wo 

to the radial velocity data. The period was fixed to 
the known optical period of 2.091 hours, and the phase 
shift 00 1 amplitude if, and the relative velocity zero- 
point vq were allowed to vary. The best-fit values are 
X = 340 ± 20 km s"! and = -1-43 ± 37 km s'^ We 
emphasize that this velocity zero point is a velocity rela- 
tive to the spectra used as our zero velocity, which were 
selected since they are the closest data to a phase of zero. 
Based on the detailed observ ation of narrow line compo- 
nent sJnJlie_2olarsJ[V Pup (jMason et al.ll2008[ ) and EF 
Eri (jHowell et al.l 120081 ) . mapping to the motion of the 
secondary star, the typical orbital phase 0.0 used for cat- 
aclysmic variables would occur near phase 0.5 as shown 
in Figures [T] and [2] However, we do not have sufficient 
data here (e.g., a velocity curve of the photosphere of the 
secondary star) to state this fact with absolute certainty. 
As these velocities are not absolute radial velocities, they 
cannot be used for cluster membership determination or 
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Table 1 

Radial velocities of EU Cnc and an unrelated neighboring G-type star 



Obs. Midpt 


Phase 


Spectrograph 


fEU Cnc,rcl 


0"(fEU Cnc.rol) 




0"(t;,,rol) 


(MJD) 




Arm 


(km s'l) 


(km s-1) 


(km s^-"-) 


(km s"-*^) 


54119.51236 


0.682 


Red 


384.9 


13.8 


27.8 


5.2 


54119.52732 


0.854 


Red 


183.0 


23.2 


9.4 


20.8 


54119.54231 


0.026 


Red 


0.0 


0.0 


0.0 


0.0 


54119.56151 


0.246 


Red 


-288.7 


27.6 


1.5 


46.9 


54119.57651 


0.418 


Red 


-187.5 


28.1 


31.5 


42.4 


54119.59146 


0.590 


Red 


201.4 


25.6 


23.8 


44.3 


54119.61542 


0.865 


Red 


179.7 


23.5 


-5.9 


71.1 


54119.51239 


0.682 


Blue 


458.5 


17.9 


36.3 


11.8 


54119.52689 


0.849 


Blue 


248.3 


27.5 


29.3 


18.8 


54119.54145 


0.016 


Blue 


0.0 


0.0 


0.0 


0.0 


54119.56153 


0.246 


Blue 


-297.7 


30.5 


9.7 


19.9 


54119.57613 


0.414 


Blue 


-158.2 


24.9 


22.9 


21.2 


54119.59067 


0.581 


Blue 


285.7 


33.6 


10.5 


22.7 


54119.61545 


0.865 


Blue 


223.2 


35.2 


22.6 


25.0 



Note. — Phase = indicates the negative zero crossing of the narrow-line component of the emission lines; velocities are relative to 
the exposure closeest to (p = 0. 



rejection. 

Qualitatively, the sine curve is not a superb fit to the 
data, especially for phases between 0.5 and 1. We note 
that, at these phases, the narrow component of the lines 
is relatively weak and the broad component strong. In 
fact, at phase — 0.59, the narrow component of the 
lines is not visible. These velocity measurements are 
likely non-Keplerian streaming motion. However, the ex- 
cellent fit for the points with phases between and 0.5, 
where the emission lines are dominated by the narrow 
component, suggests that our fit amplitude and phase 
shift are not unreasonable. 

2.2.2. Photometric Variability 

As mentioned above, absolute spectrophotometry was 
not a goal of our observations; slit losses as a function 
of wavelength could be significant and time- variable due 
to both seeing and atmospheric dispersion effects. Even 
so, it would be instructive to compare photometry de- 
rived from our spectroscopy with previous photometric 
monitoring of this system. We therefore calculate broad- 
band photometry by folding our (relative) flux-calibrated 
spectra through filter response functions. Due to the sig- 
nificant loss of UV light in the later exposures, we restrict 
our analysis to the g-band. 

In order to correct for slit losses and variable atmo- 
spheric absorption, we perform the same calculation for 
a K-type star targeted in a different slitlet on the same 
mask (note that this is not the same star used for ra- 
dial velocity comparisons in i )2.2.1[) . The star is located 
at a(J2000) 8''5r"35!46, (5(J2000) = ll°50'19'.'l, and, 
in our photometry (Williams et al, in preparation), has 
g = 20.324 ±0.032. This star also appears in SDSS DR7, 
with PSF magnitude g = 20.263 ± 0.038. Due to the 
higher precision of our data, we adopt our photometry 
for this star. 

For each exposure, the calibrated spectra of both EU 
Cnc and the comparison star are folded through the g- 
band filter response function using the Synphot synthetic 
photometry package in STSDAS. All magnitudes are cal- 
culated as AB magnitudes. Zero-point offsets in g are cal- 
culated from the comparison star; these offsets are then 
applied to the calculated magnitudes for EU Cnc. The 



Table 2 

Time Series Broad-band Photometry of EU Cnc 



Midpoint Obs. 


a,b 

y* ,mcas 


5EU Cnc, mcas^ 


SEUCnc, corr'^ 


(HJD) 


(mag) 


(mag) 


(mag) 


2454120.01789 


20.563 


20.806 


20.567 


2454120.03239 


20.676 


21.601 


21.249 


2454120.04695 


20.480 


21.672 


21.516 


2454120.06703 


20.703 


21.266 


20.887 


2454120.08163 


20.614 


20.440 


20.150 


2454120.09617 


20.467 


20.429 


20.286 


2454120.12095 


21.147 


21.706 


20.883 



^ Measured flux folded through g filter response 
^ Instrumental magnitude of the neighboring K star 
EU Cnc photometry corrected to standard system 

results are included in Table [2j The errors in this pho- 
tometry are uncertain. The random errors due to photon 
shot noise in the source and sky are small, < 0.003 mag. 
However, systematic errors such as differential slit losses 
likely dominate the photometric uncertainty. The mag- 
nitude of the error in absolute photometry is at least 
0.03 mag (the error in the broadband photometry of the 
comparison star). 

We compare our derived photomet ry with the time- 
series photometry of iNair et al.l (|2005( ). These published 
data were taken in V; we converted these magnitudes to 
g using the Populat i on I tr ansformation equations in Ta- 
ble 4 of lJordi et al* C20 06D and the B-V= 0.41 color of 
EU Cnc from [Gilliland et al.l ()1991h . As these transfor- 
mations are for single stars, and as the color of EU Cnc 
is likely changing as a function of phase, we emphasize 
that these transformations are meant to be illustrative 
only. Since the ephemeris of EU Cnc is not sufficiently 
well-determined to allow us to phase the data precisely, 
we added an arbitrary phase shift to the published pho- 
tometry. 

The results of this comparison are shown in Figure 
[21 With the exception of our observation at (j) = 0.865, 
our corrected spectrophotometry and the published time- 
series photometry agree very well in both shape and am- 
plitude. 
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Figure 3. The light curve for EU Cnc. Large filled circles are the 
spectrophotometry data from this work; error bars are not included 
but are likely at least 0.03 mag. Filled triangles with er ror bars 
(red in the online version) are data from lNair et al.l II2005I ), shifted 
in phase by an arbitrary offset. Stars (blue in the online version) 
are spectropolarimetric data form this work. The light curves are 
qualitatively nearly identical between these data sets, suggesting 
little change in the state of the polar between early 2004 and late 
2007. 



2.3. Time-series Polarimetry 

Circular polarimetry of E U Cnc wa s obtained with the 
SPOL spectropolarimeter (jSchmidt e t al. 1992) on the 
2.3-m Bok Telescope at Steward Observatory in 2007 
December. The observations were taken in the imaging 
mode of the instrument using a Hoya HA30+Y48 filter 
combination, giving a broad bandpass of ~ 4800 — 7000 
A. Image acquisition and data reduction follow those de- 
scribed in Smith et al. (2002), but modified as appro- 
priate for circular polarimetry. A A/4 wave plate is used 
to convert incident circular polarization to linear and the 
WoUaston prism separates the light into two orthogonally 
polarized beams that are focused onto a CCD. Two sep- 
arate reads of the CCD are made with the wave plate 
rotated through four positions to determine the circular 
V Stokes parameter. Two consecutive 10-minute expo- 
sures (150 s per wave plate position) were obtained on 
each of two nights. The sky was clear during the ob- 
servations, but no attempt at absolute calibration was 
made because of the non-standard filter bandpass used. 
Photometric information was extracted using a circular 
aperture of radius 6". 

The observing log and resulting circular polarization 
measurements [V/I — v) are given in Table[31 Also given 
are the background subtracted counts and the phase of 
the midpoint of each observation relative to first exposure 
on the first night assuming a photometric period of 2.091 
h. There is a significant detection of polarization in the 
first exposure, with v = —12.7 ± 1.8%. The degree of 
polarization and the total fiux both drop significantly 
in the second exposure (ti = —6.3 ± 2.3%). The third 
exposure, taken two nights later, is at nearly the same 



phase as the second exposure. The total flux is slightly 
lower and there is no significant detection of polarization 
[v = —2.8 ± 2.6); the final exposure is the faintest of the 
four and has no significant polarization. 

As these observations were taken 11 months after the 
spectroscopy, and as the ephemeris of EU Cnc is not suf- 
ficiently well known, the spectroscopic phase of these ob- 
servations cannot be calculated - the accumulated uncer- 
tainty in phase giv en the 0.002 hr uncert ainty in the pho- 
tometric period of iGilliland et alJ (|1991l ) is ~ 7.5 cycles. 
However, from the total counts in the polarimetry mea- 
surements, we know that these observations were taken 
on the declining portion of the light curve. We there- 
fore calculate relative magnitudes from the total inten- 
sity measurements in Table |3] and add an arbitrary mag- 
nitude zero point and phase shift to place these points in 
the folded light curve of Figure |3l 

From this exercise, we see that the rate of decline in the 
observed intensity ag rees with that obs erved in the time- 
series photometry of iNair et al.l ()2005f ) and in our spec- 
trophotometry, indicating that the spectroscopic phase of 
the first polarimetric observations was (/) « 0.5 to 0.65. 
Comparison with the time-series spectra in Figure [T] 
shows that this corresponds to a phase when the cy- 
clotron humps are dominant in the spectrum; these 
humps then vanish by a phase (j) = 0.85. This is con- 
sistent with the strong degree of polarization observed in 
the first exposure, and the weaker/insignificant polariza- 
tion in the other polarimetry exposures. 

In summary, we detect significant polarization from 
EU Cnc at a phase likely corresponding to strong spec- 
troscopic cyclotron emission features, and the degree of 
polarization decreased below our detection threshold at 
phases when the cyclotron features were fading or not 
evident in the spectrum. This proves that the features 
are indeed cyclotron emission. From a semantic point of 
view, this detection of significant polarization also con- 
firms that we can use the moniker "polar" to refer to EU 
Cnc. 

3. DISCUSSION 
3.1. Cluster Membership 

EU Cnc is p rojected k, I'J away from the cluster cen- 
ter defined by iMontgomerv et all (|1993D (0.4 pc at the 
cluster distance), well within the observed core radius of 
« A'.b (e.g., Tinslcy & King 1976; Bonatto & Bica 2003). 
However, this proximity alone is not sufficient evidence 
o f cluster mern b ership . 

iBellini et al.l (|2010[ ) present proper-motion data for 
white dwarfs in M67, precise to V ~ 26, obtained from 
multiple- epoch imaging from the Canada-France-Hawaii 
Telescope and the Large Binocular Telescope. Figure [4] 
shows the vector-point diagram of objects within 20' 
from the cluster center (right panel). The circle indi- 
cates the pr o per m otion membership cut adopted by 
IBellini et al.l ()2010[) : the error bars show the displace- 
ment uncertainty for EU Cnc, solidly inside the member- 
ship circle. In addition, the proper-motion selected color 
magnitude diagram (left panel) shows that EU Cnc lies 
on the WD cooling sequence of M67. 

Astrometry performed on EU Cnc and 3 encircling 
brighter known member stars, using the POSS I and II 
plates as well as a deep, '.'4 seeing image obtained by us 
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Table 3 

Polarimetric Observations of EU Cnc 



Observation Date 
(UT) 


Obs. Time 
(UT) 


Relative Phase^ 
(A*) 


Counts 
(ADU) 


VII 

(%) 


a{VII) 
(%) 


2007 December 14 


10:28 




74147 


-12.7 


1.8 


2007 December 14 


10:39 


0.08 


55485 


-6.3 


2.3 


2007 December 16 


08:41 


0.10 


48303 


-2.8 


2.6 


2007 December 16 


08:52 


0.19 


30698 


4.6 


4.0 



Phase relative to first polarimetric observation 





1 1 1 1 1 1 1 1 1 1 1^ 1. 1 1 1 1 1 1 


.1 l_ 1 1 1 1 1 1 .1 1 J 1 1 .1 -1 I- 1 ■ 1 


16 


% 




30 


V. ~ 






■ ■** ■ ■ ' \ 




32 




■ ■ ■ 1 i: ■ .;• ■ ■■■ ■Vv-Z^f '^vf-'"-'' /.- : ^; ' 


34 


. , . . 1 , . . . 1 . . . , 1 . . , . 


, ■ , ,■ ■ 1 , ,■■ , ■ !■ ,■■ -1 , . , 




10 12 3 


-0.5 0.5 1 




B-Y 


(pix/10.13 yr) 



Figure 4. The proper-motion selected color-magnitude diagram 
(CMD, left) for M67 and the vector point diagram (VPD, right) 
for all stars in the field as determined by Bellini et al. (2010). The 
cross (red in the online version) in the CMD indicates the location 
of EU Cnc. The circle (blue in the online version) in the VPD 
indicates the selection criteria Bellini et al. used to select clus- 
ter members; the error bars (red in the online version) show the 
location and uncertainty in the proper motion of EU Cnc. The 
qualitative data are convincing that EU Cnc is a proper-motion 
member of the cluster. 

at the Kitt Peak WIYN telescope in 2007, show that the 
polar's location relative to the other three stars changes 
by O'.'S ± O'.'S over the 60 year interval. This is a similar 
uncertainty to that for each of the brighter stars relative 
to each other. We therefore feel confident that EU Cnc 
is a member of the M67 star cluster. 

3.2. A low-state polar in M67 

iNair et all ()2005[ ). suggest that the luminosity of EU 
Cnc is lower than other high state polars and fully con- 
sistent with the luminosity of low-state polars in globu- 
lar clusters; the distances to field magnetic CVs are too 
uncertain to allow for such a precise comparison. The 
spectra obtained in this study, as well as those prior, 
show the clear indications of a low state polar: a steep 
Balmer decrement (compared to a flat decrement and 
Balmer jump in emission in high states) and the separa- 
ble na rrow and broad emission line components. iWarneiH 
(|1999| ) examine the high-low state range of polars and 
notes that the magnitude difference is typically near 3-4 
for a 2 hour polar. Additionally, low state polars have 
absolute magnitudes near 11-12 in V compared with 8-9 
in V when in a high state. If EU Cnc were in a high 
state with its observed apparent V magnitude, it would 
need to reside at a distance of near 3100 pc. However, 
in its low state and F ~21, we find it to be at a distance 
of approximately 850 pc or the same distance as M67 
(jSaraiedini et alJl2009D . 

Additionally, we detect the secondary star in the op- 




6000 6500 7000 

Wavelength (A) 

Figure 5. Red spectrum from phase 0.865 in Fig. 1 along with 
a scaled M5V spectrum taken from the Jacoby Atlas (star 57, 
IJaco bv et al.'lgsj'l. The M5V star is scaled 1.4 magnitudes fainter 
than the stellar continuum, approximating the level of contribution 
is provides to the red end of the spectrum of EU Cnc. 

tical spectrum of EU Cnc (see the late phase red spec- 
tra in Fig. [T] and Fig. [5]) by noting the TiO absorption 
bands chopping into the spectrum at the characteristic 
locations near 68OOA and 7IOOA. The shape and rela- 
tive amplitude of the TiO humps are a good match to a 
M5V star which is also the expected spectral type of the 
secondary star in EU Cnc (see below). Noting the dilu- 
tion of the TiO bands compared to a single M5V star, 
we determine that the secondary star contributes 20-25% 
of the continuum redward of 6OOOA. This estimate sug- 
gests that the M5V secondary star will have an apparent 
magnitude near 22.5 and, for an its My = 12.6, yields a 
distance to EU Cnc of 832 pc, again the same as M67. 
The M star continuum is also the likely cause of filling 
in the Zeeman spilt Ha absorption components. 

3.3. The masses of the white dwarf accretor and its 
progenitor star 

An advantage of finding a magnetic CV in an open star 
cluster is that numerous constraints can be placed on the 
progenitor system. As EU Cnc is a member of Messier 
67, its initial metallicity and total evolutionary age are 
identical to the cluster's characteristics. With some rea- 
sonable assumptions, we can go one step further to con- 
strain the white dwarf's mass and cooling age, which we 
can then use to constrain the white dwarf progenitor's 
mass. We describe these steps in detail below, but note 
that this methodology has been used in numerous open 
cluster and field star studies to determine WD progen- 
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Williams et al. 



itor masses (e.g. iLiebert et al.ll2005 


;lKalirai et al.l 120081: 


Catalan et al.l 120081: iWilliams et al. 


120091 iDobbie et al.l 


I2012D 



WD masses are often determined through model at- 
mosphere fits of the WD spectrum, but absorption lines 
from the white dwarf photosphere are not convincingly 
evident in our data, being filled in and misshaped by the 
emission. However, we n o te th at ILiebert fc StockmanI 
{1985) and iMason et all H2008D determined that the 
"center of mass" for the narrow emission lines originates 
between Li and the center of mass of the donor, that 
is they approximate the center of mass of the secondary 
star. We therefore assume that the narrow component of 
the emission lines originates from the center of mass of 
the secondary star and use our velocity amplitude to de- 
termine the mass function of the system and to estimate 
the mass of the white dwarf, Af-yvD- 
The mass function f{M) of the system is 

' (MwD+M2)2 2^G ^ ' 

where M2 is the mass of the donor star. Assuming a 
circular orbit (e = 0) with a period P — 2.091 hr and 
the velocity amplitude K = 340 km s^^, from Eq. 1 we 
find f{M) = 0.356Afo- 

Since we do not see ecli pses, we can constrain the in- 
clination to be < 74° fe.g. lWarneij [T995). If we assume 
MwD >> M2, we therefore find that A/wd > O.4OM0. If 
inste ad we assume an average value of < sin^ «>= 0.679 
(e.g. lTrimbl(j[T97l . then we find A%d = O.52Af0. 

Howe ver, the mas s of t he donor star is not negli- 
gible. iHowell et al.l ()200lD present detailed evolution 
models for the secondary stars in cataclysmic variables 
and note that for those systems with orbital periods 
below the period gap (< 2.5hr), the secondary stars 
follow the normal main sequence mass-radius relation- 
ship. We therefore expect the mass of the secondary 
star in EU Cue to be M2 = O.21Af0 with a radius of 
i?2 = O.22i?0. These value are roughly those of an M5V 
star. Adopting this mass for the secondary star, we find 
MwD > 0.68 Af0 (since i < 74°) and Mwd = 0.83 Mq 
for < sin^i >= 0.679. 

We note that these masses assume that our velocity 
amplitude is correct. As noted in ij2.2.1[ our value of 
K assumes that our velocity amplitude is resolving the 
narrow component of the emission lines and that the nar- 
row lines trace the center of mass of the secondary star. 
Higher spectral resolutions would be necessary to test the 
first assumption, and so we do not have good constraints 
on the errors or our confidence limits on the white dwarf 
mass. However, the WD mass estimates are similar to 
or higher than the spectroscopic masses of young mem- 
ber WDs in M67 {M = 0.5 Mq to 0.6 Mq, Williams et 
al. in preparation). Further, through a lucky coinci- 
dence, our ignorance of the WD mass does not signif- 
icantly impact estimates of the WD's progenitor mass. 
We therefore proceed (perhaps quixotically) to constrain 
the WD's progenitor mass. 

To co nstrain the progenitor m ass, we use the method- 
ology of lWilhams eFaLl (f2009l) . To summarize, we use 
the effective temperature and mass of the WD to de- 
termine its cooling age (the elapsed time since the WD 
emerged from the AGB progenitor star). This cooling 



age is subtracted from the cluster age to get the progen- 
itor star's nuclear lifetime; stellar evolutionary models 
are then used to infer the progenitor star mass. 

The WD mass estimates are sufficiently high to con- 
clude that the WD has a carbon-oxygen core, indicat- 
ing that the common envelope phase for the progenitor 
system did not occur before helium ignition in the WD 
progenitor. Since the AGB phase of stellar evolution is 
relatively short compared to the WD cooling times and 
progenitor nuclear lifetimes we estimate for this star, we 
assume that any effect of the common envelope phase 
on the calculation of the progenitor nuclear lifetime is 
minimal. 

To estimate the effective temperature of the WD, we 
assume that the minimum luminosity of the system is 
due to solely to the combined light of the M5V sec- 
ondary and the WD photosphere, and we assume that 
the WD has not been heated by the ongoing accretion. 
Under these assumptions, the WD has My — 12.5. We 
then interpolate WD photometric evolutionary tables 
provided online by P. Berge rorEI and computed from colo r 
and model calculations of 'Holberg & Bergeron' ('2006') , 
i Kowalski fc Sau mon (2006), Trcmblay et al. (2011), and 
iBergeron et al.l (|26ll[ ) to find the WD's cooling age as a 
function of mass at which its luminosity is equal to the 
observed minimum luminosity of the CV. Despite the un- 
certainty in the WD mass, the cooling age of the WD is 
relatively insensitive to its mass, at least over the range 
of O.5Af0 to I.OA//0, at logTcooi = 8.86 to 8.95. If the 
WD has experienced significant heating, or if the min- 
imum luminosity is significantly contaminated by light 
from the accretion stream and/or the accretion column, 
then its actual cooling age will be older than these cal- 
culations. This means that our derived progenitor mass 
is a lower limit on the WD's actual progenitor mass. 

We now calculate the WD progenitor mass. We adopt 
an age of M67 of 3.5 Gyr to 4.5 Gyr, and a cluster metal- 
licity of Zm67 ~ ^0 and the calculate d primordial value 
of Zq = 0.0142 (jGrevesse et al.ll2010D . Since the lower 
limit on the WD cooling age is well-constrained, we ob- 
tain a fairly precise lower limit on the progenitor mass 
of AfvvD,init = 1.43Af0 ± O.O7A'/0, with the majority of 
the uncertainty due to the uncertainty in the age of M67. 
For comparison, the current main sequence turnoff mass 
(using the same models and inputs) is 1.35Af0. 

If we relax the assumption of no heating of the WD, 
then the WD's cooling age (time since its emergence from 
the red giant progenitor) could be significantly longer; 
this would imply a larger progenitor mass. Therefore, 
we can only constrain the WD progenitor mass calculated 
above to be a lower limit; i.e. AfwD,init > 1.43Af0. 

4. CONCLUSIONS 

We have presented new observations of the polar EU 
Cue. In particular, proper motion studies strongly indi- 
cate that this system is a bona-fide member of the old 
open star cluster M67, which provides strong constraints 
on EU Cue's total age, metallicity, and progenitor mass. 
Using phase-resolved optical spectra and polarimetry, we 
have shown that EU Cue exhibits all the properties of a 
low state polar. Both assumed "best fit" polar parame- 
ters as well as standard values for the donor star lead to 

^ |http: //www. astro .umontreal . ca/- {}bergeron/CoollngModels| 
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conclusions that are consistent with the system's mem- 
bership in M67. EU Cnc is one of only four CVs con- 
firmed to reside in open clusters and the only confirmed 
magnetic CV in an open cluster. More detailed photo- 
metric, polarimetric, and spectroscopic observations on 
large telescopes will be required to refine the system 
parameters and model the system with state-of-the-art 
analyses. However models of magnetic cataclysmic vari- 
able formation and evolution are not usually attempted 
and never well constrained; further study of EU Cnc will 
greatly aid in these endeavors. 
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